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Abstract. To construct a no-arbitrage defaultable bond market, we work on 
the state price density framework. Using the heat kernel approach (HKA for 
short) with the killing of a Markov process, we construct a single defaultable 
bond market that enables an explicit expression of a defaultable bond and credit 
spread under guadratic Gaussian settings. Some simulation results show that 
the model is not only tractable but realistic. 

Keywords: (non-)systematic risk, state price density, killed HKA, Markov 
functional model, quadratic Gaussian. 



1. Introduction 

The HKA, which is an abbreviation of "Heat Kernel Approach to interest rate 
modelling, was introduced by one of the authors and his collaborators in [2J. Briefly 
speaking, HKA is a systematic method to produce a tractable interest rate model 
which is "Markov functional" in the sense of Hunt-Kennedy-Pelsser [8]. In the 
fundamental paper [2 , four different types of implementation methods are intro- 
duced. Namely, 1) Eigenfunction models 2) Weighted HKA 3) Killed HKA and 4) 
Trace Approach. As is pointed out in [5] , the eigenfunction models are tailor-made 
for swaption pricing, and a deeper understanding for its mathematical structure 
leads to the trace approach, which is mathematically most involved. The weighted 
HKA is extended to time-inhomogeneous setting and applied to information-based 
models by J. Akahori and A. Macrina pQ. 

In the present paper, we will demonstrate how the Killed HKA is applied to the 
modelling of defaultable bonds by constructing a market where the market price 
of risk and the default probability are "built in the same block" (whose precise 
meaning will be given later). We stress that the HKA is basically a state-price 
density approach where everything is written under the physical^ statistical mea- 
sure. Since the HKA furthermore gives an analytically tractable model in nature, 
the framework proposed in this paper would be promising in respect of modelling 
defaultable markets. 

The organization of the present paper is as follows. After recalling the plain- 
vanilla HKA in section [2~T1 and the killed HKA in section [2~2l we shall give the main 
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result, a framework with in the Killed HKA to model a defaultable bond market in 
section |31 In sectional we will give some simulation results of an explicit example 
with a quadratic form of Wiener process. 
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2. Heat Kernel Approach 

Here we briefly recall the approach. 

2.1. Plain- Vanilla HKA. We work on a probability space (f2, J 7 , P) with filtration 
{J-"t} t>0 . Now we consider a general Markov process {Xf } t>0 xeS on a polish space 
S. ~ " ' 

Definition 2.1. Let X be an S -valued Markov process. We shall say that a function 
p satisfies the propagation property if 

(1) E\p(t,X*)]=p(t + s,x) 
holds for any t, s > and i£5. 

The following fact is initialized by [3] and developed in [2]. 

Proposition 2.1 (Akahori et al. [2]). Let X be a S-valued Markov process, X be a 
positive function on the half line, and p be a function with the propagation property. 
The bonds market given by 

(2) P(tT)- P( X T + T-t,X?) 
is an arbitrage free market. 

Example 2.1 (Generic Example). Take a measurable, bounded h : S R>o. then 

(3) p(t,x):=E[h(Xf)] 

satisfies the propagation property {!]]. In fact, by the Markov property, we have 

E [p{t,X s )\ = E [E [h{xf>)]] = E [E [h{Xf +s )\T^]] 

= E[h(Xf +s )] =p(t + s,x). 

Proof. Let 7r 4 = p(\ t ,Xf). By the propagation property of p and the Markov 
property of X, we have 

E[n T \T t ] _ E[p(A T ,Af)|J- t ] 



Pf(t,T) 



E 



Tit P(^t,X. 

'p(\ T ,X*l t ) 



p{\ T +T-t,Xf) 



p{X u Xf) p(X t ,Xf) 
This means tt is the state price density of the market. □ 

It should be noted that we do not assume Tr t to be a supermartingale in |2], i.e. 
in economic terms we do not assume positive short rates. The four implementation 
methods mentioned in the introduction is introduced in [2] to obtain supermartin- 
gales out of a propagator, or equivalently to obtain positive rate models. 



DEFAULTABLE BONDS VIA HKA 



3 



2.2. The Killed HKA. We then recall, and give a more detailed description to, 
the Killed HKAQ. Let V be a non-negative measurable function on S. Put = 
y + J Q V(Xg) ds for y e R. Let us define 

(4) <j(t,x)=E[exp(- / ^(X s K )rfs)]. 

Jo 

Then the function 

q(t,x,y) = e~ v q(t,x), 
satisfies propagation property with respect to (X X ,Y V ); 

(5) E[q(8,Xt,Y?)] = q(t + 8,x t y). 
In fact, by the Markov property of X, we have 

E[e- yt + s |J-t] = E[e~ YA6t) \T t } X e~ y ' = E[ e - r '|X t ] x e^, 
where 9 is the shift operator. Thus we obtain 

E[q(s,X?,Y?)} =ne- Y * V q(s,X?)} = E[ e - y *VW] = g (£ + S , x, y). 
This fact ensures that the bond market model constructed as 

feA ph t\ g(Ar + r-*,^f) 

(6) ^' T)= ? (A t) X f ) ' 

where A is an increasing function, is arbitrage-free since we can choose 



n t =q(\ t ,X t )exp(- f V(X s )d S ) 
Jo 

as a state price density of the market. In fact, 

E[7T T | -Ft] = E[ E[exp (- f + T V(X£) du)\T T ] exp(- / F(X S ) ds) |JT t ] 

Jo 

= E[exp(- f + 1 V{XZ)du) |X t ]exp(- f V(X s )ds) 
Jt Jo 



'/(A/ +7 -/.A'/l.-x].!- / V(X s )ds)=ir t ^' +T LXi) 



q{Xt,X t ) 

Note that the bond price P is decreasing in T since q is increasing in £, which is 
ensured by the positivity of V. Thus we obtain a positive rate model. 



3. HKA to Defaultable Bond 

This section is the main part of the present paper. Let us now consider a de- 
faultable bond in the following situation: 

(1) The bond pays a unit account at the maturity T unless it defaults. 

(2) At the default time r, nothing will be recovered. 

(3) The state variable is a Markov process {Xf;t > 0}, which is observable in 
the market. 



^This part is shared with [9]. 
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(4) The default probability is completely determined through the information 
of X in the following manner; the hazard rate of the default time on the 

filtration T x is given by E[l{ r>t} | = cxp (- J V(X*) du^j where T x 

is the natural filtration on X and V is a non-negative measurable function. 

(5) The default come as a "surprise" to the market. To be precise, the market 
nitration {Q t } is denned as Qt = a(X s ,{r < s};s < t) and assume that 

/ /-At 

(6) A state price density of the market is given by 7Tt := q(\t,X t ) = E 



expl -J V(JT*)du 



where q is defined as @ and A is a non-decreasing function. 

Note that the assumptions 1-5 may be natural (except assumption 2, which 
assumes zero recovery) and very generic, while the last assumption is very specific 
in that the function V controls both the market price of a risk as well as the default 
probability of a bond. Very heuristically speaking, this market is fully subject to 
the risk of a defaultable bond. 

We stress that this is just a toy model, which exhibits how the killed HKA is 
applied to a defaultable market modeling. The following is established in [9]: 

Theorem 3.1. Under the above assumptions 1-6, 

(i) the price Pd(t,T) of a defaultable zero coupon bond is given by 

(7) P d (t,T) - l {T>f} ^— ^ , 

(ii) the price Pf(t,T) of a default-free bond is given by 

q(\ T + T-t,T-t,X?) 



(8) Pf(t,T) = 
where 

(9) q{t,s,x)=E 



g(At,Xf) 



cxp (-J V{XZ)du 



(iii) and then the "credit spread" is given by 
(10) ' Tl0S q(X T + T-t,0,Xn ■ 

Remark 3.1. Note that the "credit spread" makes no sense when r < t, so we can 
only think of the case that t > t. 

Proof. The proof is based on the following fundamental lemma due to Dellacherie 
(see [6]): For any J^-integrable random variable Z and < t < T, we have 

Hr>t} 



ni {T>T} z\g t ] = ^ E[i {T>T} z\m 

^i 1 {r>t}Kt I 



Hence, we have 



P d (t,T) = 1 1{T>t} E [l {T>T} n T \^] . 

n t & [ l {r>t}\^t J 

Then by a Markov property and a Tower property, 

E [l {T>TFT |jf] = E [ E [1 {T>T} |^] q(\T,XZ)\r?] 
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and by the assumption that E[ l{ T >t} I F*] — ex P ^~ J V(X^) du 
E [ E [ 1 {T>T} | F$] q(X T ,XZ) | jf ] = exp (- [ V(X*) du) E \e~ F V W> ds q(\ T , X$) \ T? 



Here applying the fact of the equation ((3J), we have 



E 



so that 



e- v ^ ds q(\ T , X*)\F t x ] = g (A T + T — t, Xf), 
P(tT) i ?(AT + T-t,*f) 

F d {t,l ) = l{r>t}" 



3(At,Xf) 

On the other hand, (ii) follows a Markov property and a Tower property. And it is 
known that the "credit spread" is given by 

-or log ^> 
g P f (t,T) 

Here since Pd(t,T), Pf(t,T) are (i), (ii) respectively, 



-dr log 



Pd(t,T) 



dr log 



E[g(A T ,Xf)|^ 



vi 



P/(t,T) °g(A T +T-t,Xf) 
when t > t. Then by a Markov property and a Tower property, 

E[g(A T ,X|)|jf] = E 



Hence, we obtain 



-Or log ; — -rr = Or log 

& Pf(t,T) s 



:g(A T + T-t,T-t,Xf). 

g(A T + r-t,r-t,xf) 

? (A T + T-t,Xf) 



□ 



4. Quadratic Example 
Now we give some simulation results of an explicit example, where X is a d- 

o2 |2 

dimensional Wiener process and V(x) = s -f L - (/3 > 0). Let q(t,x) and be 
as in (dJ), (O, then they are explicitly given by 

pi, ?(t , lH(OThW -^ P (-^^), 

and 
(12) 

«(,,,) . ( - .) + ^»„h« t - .) exp (-^ ^X^W-J ' 

which result from Lemma |6.1[ Corollary 1 6 . 1 1 in the Appendix. Hence, we obtain the 
analytic expression of the bond prices. The following simulated yield curves (FigfTJ) 
and (Fig|2Jl implied by a default-free bond as 

- i log P f (t,T) 

are obtained by using © and (|12p. Here the parameters are set to be (3 = 0.1, 
x = 0.01, 10, 20, 30, A t = e*/10, and the present time t = in (FigHJ), /3 = 1.8 and 
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Figure 1. Simulated 
yield curves implied by a 
default-free bond when 
At = e'/lO 



Figure 2. Simulated 
yield curves implied by a 
default-free bond when 
A f = e'/lOO 



(FigfT]) and (Figj2j) show increasing the value of x does not make the curve shift 
upward, but also cause a "hump" in the curve, which can not be observed in the 
normal afhne model, with the proper choice of A. Moreover, using the formula (|10j) 
and ((HI), we obtain the following simulated credit spread curves as (FigJ3]). Here 
the parameters are set to be At = Vt, P — 0.1, 0.2, 1, x — 0, and the present 
time t = 0. 



Figure 3. Simulated credit spread curves 



As usual, the lower the credit rating of a defaultable bond is, the wider the 
spread is, and it is non-decreasing in the maturity time. Moreover, the spread of a 
defaultable bond lower rated is much wider in the maturity time than the one of a 
bond higher rated. It should be thought that (FigEJ) shows this fact. 
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5. Conclusions 



We have introduced a way of constructing a single defaultable bond market model 
under the physical measure P by applying the killed HKA. We have also presented 
some simulation results in the quadratic case. Comparing the well-known Hull- 
White model, we can have observed a complex "hump" in the yield implied by a 
default-free bond, which comes from the parameter A. 

6. Appendix 

Lemma 6.1. Let X be a d- dimensional Wiener process starting at x. For a,0>O, 
it holds 



(13) 



E [ e -a|X?| a -^/*|x;| a d.] 

'cosh (3t + ?f sinh fit 



-d/2 



exp 



(2ai + l)- d / 2 exp(- 5 H f d_ 



(-¥ 



sinh (3t — 2a cosh (3t 
(3 cosh fit+2a sinh fit 



p>o, 

13 = 0. 



This is well-known formula and there are many ways to prove it. One way is 
presented in [5]. 

The following is an immediate consequence of Lemma 16.11 



Corollary 6.1. Let X be a d- dimensional Wiener process starting at x. For /3 > 0, 
it holds 



E 



= (cosh P(t-s) + (3s sinh /3{t - s)) d/2 
Px 2 sinh/3(t-.s) 



x exp 



2 cosh/3(t- s) + /3ssinh,S(i- s) 
Proof. By a Markov property, a Tower property, and Lemma 16. 11 



E 



E 



E 



■/*|XJ| 3 d«|j7X 



E 



E 



cosh /3(t — s) 



-d/2 



E 



exp 



/3|Xf| 2 sinh^(t- s) 
2 cosh/3(t - s) 



The proof is complete by replacing a by ^ 



sinh P(t-s) 



2 cosh^(t-s) 
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